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We derive the one loop renormalization group equations for the Cabibbo-Kobayashi-Maskawa 
matrix for the Standard Model, its two Ifiggs extension and the minimal supersymmetric extension 
in a novel way. The derived equations depend only on a subset of the model parameters of the 
renormalization group equations for the quark Yukawa couplings so the CKM matrix evolution 
cannot fully test the renormalization group evolution of the quark Yukawa couplings. From the 
derived equations we obtain the invariant of the renormalization group evolution for three models 
which is the angle ct of the unitarity triangle. For the special case of the Standard Model and 
its extensions with v\ « V2 we demonstrate that also the shape of the unitarity triangle and the 
' Buras-Wolfenstein parameters p — {l — |A^) p and f) = {l — ^A^) are conserved. The invariance 

of the angles of the unitarity triangle means that it is not possible to find a model in which the 
■ CKM matrix might have a simple, special form at asymptotic energies. 
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I. INTRODUCTION 



^ To obtain finite results in quantum field theory, in a higher order than the tree level, one has to perform the 

renormalization program. The independence of the renormalization procedure of the renormalization point leads to 
' the dependence of the Lagrangian parameters on the point of renormalization. This dependence is governed by the 
renormalization group equations for the coupling constants and other parameters of the Lagrangian [ll,!^. The most 
<^ I important application of the renormalization group equations is a possibility to study the asymptotic properties of 
the Lagrangian parameters like the running masses and the coupling constants (for a very early example see e.g., 0]). 
^\ , The theories at the asymptotic energies may reveal some new symmetries or other interesting properties that give a 

■ deeper insight into the physical content. Also from the physical requirement of the stability of the theory one can 
e.g., estimate the range of the physical parameters like the Higgs mass in the Standard Model or its extensions [J| 

■ and the references therein. 

T— I ' The observable parameters of the Standard Model are the following: 6 quark masses, 3 lepton masses, 4 parameters 
OO . of the Cabibbo-Kobayashi-Maskawa (CKM) matrix Q and 3 gauge coupling constants. The CKM matrix is 

' obtained from the diagonalizing matrices of the Yukawa couplings. 
^ , The renormalization group (RG) equations for the CKM matrix is obtained from the renormalization group equa- 
tions for the Yukawa couplings. In this paper we derive the RG equations for the CKM matrix for the Standard Model, 
. , two Higgs doublet extension of the Standard Model and the Minimal Supersymmetric Standard Model. We first get 
' a complete one loop equations without any approximation and next discuss these equations taking into account the 
hierarchy contained in the Yukawa couplings and the CKM matrix elements. 

The interesting fact that can be observed is that the equations for the CKM matrix do not depend on all the 
parameters of the original equations for the Yukawa couplings. Finally we discuss the physical implications of the 
derived equations. 
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II. RG EQUATIONS FOR THE CKM MATRIX 
A. Preliminary considerations 

The part of the Lagrangian of the Standard Model and its extensions that contains the Yukawa interactions of the 
quarks has the foUowing structure 



(1) 



Here the jju and the yd are 3x3 matrices of the Yukawa couphng of the up and down quarks, $1 and $2 are the 
two Higgs doublets (in case of the Standard Model $2 = $, <&i = e$*, e is the 2x2 antisymmetric tensor). The 
qL and {qu,d)R are the left-handed quark doublet and right quark singlets, respectively. Eq. ([1]) implicitly contains 
the summation of the quark generations. The Yukawa couplings can be diagonalized with the help of the biunitary 
transformation: 



•lu - uIy^Ur, yd - D{YdDR 



(2) 



where Ul, Ur, Dl, Dr are unitary, diagonalizing matrices and Yu and Yd are the diagonal matrices with positive 
diagonal elements. The quark masses (running) are equal to 
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^Wl(l^n)ii, mc = -j=Vi{Yu)22, 
-^W2(>d)ii, m, = -^W2 (>d)22 : 
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1 , , 

{Yd) 33 . 



(3) 



Here Vi and V2 are the vacuum expectation values of the Higgs fields and in the Standard Model vi ~ V2 — v. 

From Eqs. ^ one can see that there is a hierarchy between the elements of {Yu)u and bewteen the elements of 
{Yd)ii that follows from the hierarchy of the quark masses. 

The CKM matrix of the charged weak current is obtained from the left diagonalizing matrices Ul and 



(4) 



The matrices Ur and Dr are not related to any observables in absence of the right handed currents. 

The general structure of the one loop renormalization group equations for the Standard Model and its extensions 
for the Yukawa couplings yu and yd is the following 0, [1] 



dyu 
dt 

dyd 
dt 



1 



(47r)2 
1 



i^tit) + ^yuyl + «3Tr(y„yI) + ^iivdyl) + <Mydv\)]yd- 



(5a) 
(5b) 



Here t = ln(/i/mj), is the energy of the renormalization point and nit is the top quark mass. Eqs. ^ form a system 
of coupled non-linear equations. The parameters Q;"(t), ctfit) are quadratic functions of the gauge coupling constants 
and , af , i = 2, . . . , 5 are constants. The values of the af{t) and af{t) do depend on the model and they are given 
in Appendix. 

B. Derivation of the RG equations for the CKM matrix 

We will now obtain from Eqs. ([5]) the renormalization group equations for the elements of the CKM matrix 
The derivation is done in the following steps: 

1. Insert representation Q in Eqs. ([5]). 



^ We use the following notation 



V21, V22, V23 
V3I, V32, V33 



Kid: Vui 
Vtd, Vu 
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2. Eliminate the matrices Ur and Dr from the equations. 

3. Derive the final form of the equations for the CKM matrix. 

Before the derivation, let us observe that after the differentiation of the unitarity relations U^U = 1 and JJJ/^ = 1 we 
obtain 



dU 



dW fdU 



dt 



dt 



— [/^ + U— = —U^ + — c/t 



dt 



dU, 



dt 



U 



dt J 



dt 



(6a) 
(6b) 



i.e., the matrices ^^U, and U^j^ are antihermitian. 



In the derivation of the RG equations for the CKM matrix we will explicitly show all the algebraic operations that 
are necessary for Eq. ([5a|) . the analogous steps for Eq. (j5bp are not shown and only the final results are quoted. 
Step 1 

From Eq. ^ and ([5a|) we obtain 



(«iW+«2t^i>;'t^L + "3Tr(i;2)+<( 4y/Di) + a^Tr(y,2)j uIy^Ur. 



dt (47r) 
After differentiation, the left hand side of Eq. ([7]) becomes 



dU} fdYu + dUn 

dt ^ dt ^ dt 



(7) 



(8) 



Step 2 

Using relation (O we multiply Eq. ^ from the left by Ul and from the right by Ul/Yu 



dU 



Lv 2 



dt 



dY, 
dt 



■Yu + r„ 



dU 



«rrt 



dt 



U'rYu - 



(47r) 



+ o^Wu ' + «3Tr(r„ ^) + al VY, + a^Tr(y, ^)) Y^^ (9) 



and next we add to Eq. ([9]) its hermitian conjugate. Then from the property (|6bp the term with the Ur matrix cancels 
out and we obtain 



dY. 



dt 



dul 

dt 



^ (2 {al{t) + ^ ««Tr(y„2) + a«Tr(r/)) Y^^ + al {VY^^V\Y^^)) . (10a) 



Here [ • , • ] and { • , • } are the commutator and anticommutator of two matrices. From Eq. (jSbp we get 

1 



dt 



dD\ ^ 



(47r) 



(2 («f (t) + a;^Tr(i; 2) + 2 ^ ^.Tr(y, ^)) + {V^Y^ V, }) . (10b) 



Eqs. (jlOp are matrix type differential equations of dimension 3x3. The equations for the diagonal elements give the 
renormalization group equations for the diagonal elements (Fu)^^ and (Fd)jj • We will not consider these equations 
here. The non-diagonal Eqs. pop can be written in the following form 



(4^) 



2 '-^4 -'-'ij 



D, 



dDl 
' dt 



(11a) 
(lib) 



with no summation over The matrices i?"''^ are the right hand sides of Eqs. (|11|1 and i?,"''^ are equal to 



(12) 
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The right hand side of Eqs. pT|) is antihermitian, this guarantees the unitary evolution of the matrices Ul and Dl. 
Step 3 

Now we derive the renormahzation group equations for the elements of the CKM matrix V . 
With Eqs. HI]) we obtain 



Dl = -R"--V, Ul—7^ = V-R''. (13) 



After taking the sum of two equations in (|13p we obtain 

dy = dU,^l + U,'^A. = -R--V + V-R^^ (14) 
dt dt ^ dt ^ ' 



or written explicitly in terms matrix elements with indices 

fdV 



- 



V dt 



Eqs. HH) and (fT5|) are the RG equations for the CKM matrix, but they are not yet in the final form, because as seen 
from Eq. (fTS)) they contain the diagonal elements of the matrix (i?"''*)^^, which cannot be derived from Eqs. (|lip . 
they are unknown and Eq. (jl4p cannot be directly solved. It turns out, however, that the problem of unwanted terms 
can be overcome for the absolute values of the matrix elements of the CKM matrix which we will obtain in the later 
part of the paper for the approximate equations. The exact form of these equations is not very useful. 

III. APPROXIMATIONS 

A. Hierarchy of the Yukawa couplings 

Let us notice now that Eqs. ((H]) and ([15]) are one loop equations i.e., the two loop and higher order loop terms 
are neglected. The loop order parameter of the RG equations is l/(47r)^ ~ A'^, A ^ \Vi2\, so it means that we have to 
keep only the terms of the relative order lower than A^ in Eqs. and (fT5)) . We have the following hierarchy of the 
Yukawa couplings that follows from Eq. ((3]) and the experimental values of the quark masses [§] 



Using the one loop approximation and Eq. (|16p we have 



(Y 2) 




(^u^)33 






O^d ^)33 





(16) 




H:^-'^{+1 ^or^<J (17) 



and Eq. (|15p can be written in the following form 
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1 


dVn 


Vn 


dt 


1 


dV33 


V33 


dt 


1 


dVi3 


Vi3 


dt 


1 


dVsi 



« - Rii) - {{YDii - iVYiV^),r) - j^ai {{Y^U (Wy„V)n) , (18a) 

(i?33 - ^33) + {iYi)33 - (VYiV^h,) + jl^ai {{Y^U - (V^Y^Vh,) , (18b) 

(i?33 - ^11) - {{YlU - (VYd'V^U) + j^ai {{Y^U (V^Y^VU) , (18c) 

« - ^33) + 7A2 < {{Yihi - {VYiV%,) - -^ai ((Y„2)33 _ iV^Y^VU) ■ (18d) 



V31 dt ' "33. ' (4^)2-4 v.- a/.. ^ — a■ /oo; ^4^^ 



^ We explicitly write the RG equations for Vii, ^13, V31 and V33, because these equations are sufficient for the determination of the 
evolution of all the observables of the CKM matrix, which are determined by 4 parameters. The equations for the 5 remaining elements 
can also be derived, but they are more complicated and we do not consider them here. 
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Constants of the RG evolution 



Let us notice that in Eqs. (|18[) on the right hand side the terms i?,^''' are purely imaginary and the rest of the 
terms are real. From this we can obtain the quantities that do not evolve with energy. Observing that the following 
expressions are purely real 



^ - In , (19a) 

Vn dt Vl^33 dt ) V^i3 A ) V31 dt dt V 1^31 1^3*3. 

1 dVi, ( 1 dV^^X 1 ^^^13 , ( 1 dV:,iX _ ^ lnf ^"^3V | ^^g^) 



Vii dt V^33 dt J Vi3 dt V^3i dt J dt V ^13^33 
1 dVii 1 dV:-53 / 1 dV"i3\* / 1 dV:ii\* d 



1/ ■ T/ ru ■ \M I / ,111(^11^33^^1*3^^3*1), (19c) 

Kii dt 1/33 dt \ Viz dt J \ Vzi dt J dt 

we get from Eqs. (fro| that the following functions of the CKM matrix elements are constant during the RG evolution 



r j j ^ r J J = (V^nl^33Fi*3^3*i)) = const. (20) 

We thus see that though Eqs. cannot be solved they give us an information that certain functions of the CKM 
matrix elements do not evolve. 

The first constant is equal to the angle a of the unitary triangle, the second one is equal to one of the angles of 
the unitary triangle that is obtained from the multiplication of the first and third row of the CKM matrix. The 
third constant is the phase of the rephasing invariant {VuVszVi^jV^i) , whose imaginary part Im(Vii V33V]*3V3\) is the 
Jarlskog invariant [lO|. The physical meaning of this constant is discussed in the next section. 



C. RG equations for the moduli of the CKM matrix elements 

Let us now find the RG equations for the squares of the absolute values of the CKM matrix elements |ViiP, IV33P, 
|Vi3p and [Vsip. One can easily show that the following relation holds 

^«-2Rer±^V (21) 



\V^j\^ dt \V,j dt 



Now, using the fact that R"^'^ are purely imaginary we obtain from Eqs. (|18p 

\vl/^Y - - (i^< ((^'^')" - (^^-i^^^)") - (i^"' ((^« )" - (y'^uVU) , (22a) 

^ ^ < ((1^.^)33 - (Vy,Vt)33) + -^a^ ((1^)33 - (Wr,V)33) , (22b) 



1^3312 dt (47r)2 4 ^•^■^ ^ d '■'■^^ (47r 

1 ^1^131' 2 ,....2,.t^ N . 2 



< {iYl)33 - {VY,'V^)n) + j-^at {{Y^)n ' (V^Y^VU) , (22c) 



r^^^ - (4^< - (V^>^i^^^)33) - ((i^„^)33 - {y'YlV)..) . (22d) 

These equations form the complete set of the RG equations for the CKM matrix that do not contain unnknown 
functions. Eqs. ((22)l can be rewritten in terms of the squares of the absolute values of the CKM matrix 



1 dW^if 2 



^< (|Vi3|' {Y,f + 1^12^ {Ysf - (l - \V^if) {Y,f) 



dt (47r) 

^ (47r) 



+ vA^al (IV31I' {Ytf + \V2i\' (Y^f - (1 - iFiil') {Y^f) , (23a) 



1 dlVn f 2 

|2 



1^13 r dt (47r) 



— < (- (1 - 1^131') [Y.f + \V,,\' {Y,f + \Vnf (Y.f) 

2 

(4^ 



+ j^ai IF33I' {Y,f - 1^231' {Ycf + (1 - 1^131') {Yuf) , (23b) 
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1 d 1^31 1' 2 



- (^"2 ((l - l^3i|') (Y^f - 1^211' (Fc)' - lliil' ,) (23c) 



1 clW..? 2 



^33 1' (4 



^< ((l - 1^331') {Y,f - 1^321' (K,)' - (i^.)') 



+ 74^«2 ((l - 1^331') {Ytf - 1^231' {Ycf - |V^13|' (l^.)') , (23d) 



where we use an intuitive notation {Y^) = {Y^)^^ - Y^ , {Y^)^^ = Y^^ , (r|)^^ = (rl)^^ = Y^, (Y^)^^ = 
Y^ . We can see from these equations that one needs to know the running of the [Y'^ A to obtain the evolution of 
the CKM matrix. 



IV. RG EQUATIONS FOR VARIOUS MODELS 

The hierarchyof the Yukawa couphngs, Eq. ^ and the hierarchy of the CKM matrix, best seen in the Wolfcnstein 
representation [11] permit to obtain the approximate, simpler form of Eqs. (|23p . that is compatible with the one loop 
approximation. 

For the approximate form of the renormalization group equations we will consider two scenarios: 

1. Standard Model or two Higgs extensions of the Standard Model and Minimal Supersymmctric Standard Model 
with vi ~ V2. 

2. Two Higgs extension of the Standard Model and Minimal Supersymmctric Standard model with ^ = 

The approximation principle that we will follow consists in keeping in the RG equations the terms of relative order 
order of A'^ or lower, what is also consistent with the Wolfenstein parametrization of the CKM matrix. 

In the first scenario the squares of the Yukawa couplings of the down quarks can be neglected in comparison to 
those of the up quarks. Also the squares of the Yukawa couplings of the charm and up quarks can be neglected in 
comparison to the top quark Yukawa coupling and thus the resulting renormalization group equations are considerably 
simplified and read 

' %c.i\V,,\'Y,\ (24a) 



iV^ii 


|2 dt 


1 


d\Vi3\^ 


1^13 


|2 dt 


1 


d\V3if 


\V31 


|2 dt 


1 


d 1^331' 



(4^) 



(47r 
2 



^ «2 1^^331' i;', (24b) 



l2' 



(47r)2 
2 



l-l^3l^r,^ (24c) 



These equations were derived earlier and were explicitly solved in Rcf. [H [H] so we will not discuss them in more 
detail. 

The second scenario is physically very interesting if we assume that the origin of the hierarchy of masses between 
the top and bottom quarks lies in the hierarchy of the vacuum expectation values of the Higgs fields: 

vi rrit , , 

— « — (25) 
V2 nib 

from which we obtain 

(K)33 « {Yd)^^ , (26) 



7 



so the Yukawa couplings of the down quarks have to be kept in the equations. From the hierarchy in Eq. ITBl) one 
derives after straightforward calculations the following equations for the CKM matrix elements 

1 d\V^^f _ 2 f,. 



Fill' dt (47r) 



^ («2 \V3i\' + al 1^131' + «4 1^421' , (27a) 



1 dl^isl' 2 



1^13 1' (4^) 



— (ai 1^331' + (l - 11^131') n') , (27b) 



' ^("^(l-|^3l|>*^ + <|1^33rn^), (27c) 

(a^y,2^<n2) (1-1^331'). (27d) 



1/31 1 2 dt (4 



1 d|V^33l' 2 



1 1/33 1 2 (4^^) 

Eqs. ((24| form a complete set of equations for the CKM evolution for the first scenario and Eqs. ((27)) describe the 
CKM matrix evolution in the second scenario. Eqs. require the knowledge of the evolution of the top quark 
Yukawa coupling Yt(t) while for the second scenario we need the knowledge of the Yukawa coupling evolution of the 
top, bottom and strange quarks, Yt{t), Yb{t), Ys{t). 

A. Solution for the first scenario 

Eqs. p4)) have been analized before [l^l but we shall make some additional comments here. Let us notice that after 
subtracting Eq. ([Mc]) from Eq. ((Ma|) and Eq. ((24bl) from Eq. ((24d)) we obtain 

1 d|Fn|' 1 d|y3i|'_ 2 ,^ 



o ^ 7 ^ = T:r^aiYt^ (28a) 

1^11 1' dt 1^31 12 dt (4^)2 2 * y I 



1 d 1^331^ 1 d 1^131^ _ z ^,^2 

2 J+ iT^ |2 



a^Yt^ (28b) 



1^331 '^^ 1^^131 '^^ (4^)' 

Viil^ 1 IVaal^ 

same. Moreover, if we substract Eq. (j28bp from Eq. (|28ap we get 



Eqs. (12811 can be easily integrated and from them it follows that the evolution of the ratios ' / .a and ' ^^U is the 

^ o I V31 1 I via 



1 1 d|T/3l|' 1 rf|l/33|' , 1 rf|Vl3|' 



2 |T^j2 ^-j- ,y J2 ' ,2 ^^ (29) 



IViiT 11^31 r l'^33r 1^131 

SO in the case of the first scenario the ratio 

11:^1!:^= const (30) 

|^33n^3l|' 

is also constant. 

Next from Eqs. (l24a|) . (l24cll and Eqs. (l24bl) . (pid)) one can derive that the following ratios are constants during the 
RG evolution 

1^131' , l^lll' , 

2 — const, 2 = const. (31) 



1-11^331 ^-\V3 



Eqs. (|28p and ([31]) allow to derive the exact evolution of the CKM matrix for the first scenario. From these equations 
it also follows that the evolution of the CKM matrix depends only on one function of energy. 



Discussion of the second scenario 



The analysis of the RG equations for the second scenario Eqs. ([27)1 is more involved. One can see that Eq. (|27dp . 
can be rewritten in the following way 



d /, 11/33 1' \ 2 
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so the evolution of IV33I is explicitly known. 

Next, using the solution of Eq. (j32p one can solve Eqs. (l27bll and for \Vi^\^ and |V3i|^ Eq. can be 

solved using the solutions of Eqs. (|27bp . (|27cP and (j27dp . The full phenomenological analysis of Eqs. (p7)l is beyond 
the scope of this paper and will be published elsewhere. 



V. GENERAL DISCUSSION AND CONCLUSIONS 



We have obtained the full set of the one loop RG equations for the CKM matrix for three models: Standard Model, 
two Higgs extension of the Standard Model and the Minimal Supersymmetric Standard Model. They are given in 
Eqs. ((22)) or ([23)1 and they form a set of coupled, non linear ordinary differential equations for the squares of the 
moduli of the CKM matrix elements. Eqs. (HH) were derived from the RG equations for the quark Yukawa couplings 
by the elimination of the right diagonalizing matrices. One can see that Eqs. ()22|) depend only on a\ and i.e., only 
on a subset of a"' , i = 1, . . . , 5 of the model dependent parameters of the RG equations for the Yukawa couplings, 
Eqs. ([5]). This means that the evolution of the CKM matrix can test only partially the evolution of the Yukawa 
couplings. 

The most important result of the paper is the demonstration from Eq. (|20p that the angle a of the unitarity triangle 
is constant during the evolution. The angle a is also equal to the one of the angles of the unitarity triangle obtained 
by the multiplication of the first and third rows of the CKM matrix and is also equal to the phase of the rephasing 
invariant ViiV33T/]*3^*ii whose imaginary part is equal to the Jarlskog invariant of the CP violation. For the general 
case of the three models this is the only constant that can be obtained. On the other hand for the first scenario, i.e., 
for the Standard Model and its extensions with vi « we additionaly have the constant in Eq. (|30p . which is equal 
to the ratio of two sides of the unitarity triangle that are adjacent to the angle a. From this it follows that during 
the RG evolution the unitarity triangle is only rescaled and its shape does not change. This means that for the first 
scenario it is not possible to construct an asymptotic model with some simple, special form of the CKM matrix. 

Let us observe that the constants in Eqs. ([20]) and ([30]) can be written in terms of the modified Buras-Wolfenstein 
parameters [3] P = (1 ^ l^)p a^nd ^ = (1 — \^ /2)r] in the following way 



Im In 



13 



■n 



V 



p-(l-iA2)(p2+^2) p(l_p)_j^2 



= const. 



(33) 



and 



^11^*3 



V31V, 



33 



(l-p)2 + ^2 



const. 



(34) 



Eq. ((33|) is valid for both scenarios, so we see that a relatively complicated function of the Wolfenstin parameters is 
a constant of the RG evolution for the general case. Eq. ([M]) is valid for the first scenario so from Eqs. ([55)1 and 
it follows that for the first scenario the parameters p and fj are constants of the RG evolution. 
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APPENDIX: PARAMETERS FOR THE VARIOUS MODELS 



TABLE I: Coefficients a[, l = u,d for various models. 





SM and DHM 


MSSM 


aUt) = 




, 13 2 2 16 2n 
-(Y^5i +3^2 + yPa) 




/ 1 2 ,9 2 , o 2x 

-(^ffi + -^92 +833) 


/ 7 2 , 2 ,16 2n 
+352 + Yffs) 



TABLE IL Coefficients fc = 2, . . . , 5, / = m, d for various models. 



I 


I 

<y-2 




ai 




u 




3 




3a 


d 




3o 




3 



(a, 6, c)sM = (1,1,-1) 

(a, 6, c) DHM = (0, 1, ^) 
2 

(a, 6, c) MSSM = (0,2, -) 
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